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Abstract: We discuss the renormalisation properties of the complete set of AB = 2 fom'- 
quark operators with the heavy quark treated in the static approximation. We elucidate 
the role of heavy quark symmetry and other symmetry transformations in constraining 
their mixing under renormalisation. By employing the Schrodinger functional, a set of 
non-perturbative renormalisation conditions can be defined in terms of suitable correlation 
functions. As a first step in a fully non-perturbative determination of the scale-dependent 
renormalisation factors, we evaluate these conditions in lattice perturbation theory at one 
loop. Thereby we verify the expected mixing patterns and determine the anomalous dimen- 
sions of the operators at NLO in the Schrodinger functional scheme. Finally, by employing 
twisted-mass QCD it is shown how finite subtractions arising from explicit chiral symmetry 
breaking can be avoided completely. 
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1. Introduction 

The oscillations among pairs of neutral S-mesons provide crucial information for pinning 
down the elements of the Cabibbo-Kobayashi-Maskawa (CKM) Matrix that are associated 
with the top quark. Owing to the absence of flavour-changing neutral currents in the 
Standard Model, these oscillations are described by box diagrams, in which the flavour 
change is brought about through the intermediate propagation of a virtual top quark. 
By integrating out the VF-boson, the box diagram is replaced by an effective point-like 
interaction vertex associated with the left-handed AB = 2 four-quark operator 

Oll(x) = 6(x)7^^,(x) b{x)^f:Mx) , (1-1) 

where 7^" = 7^(1 — 75), and the flavour label £ denotes either a (i or an s quark. The 
matrix elements of Oll between S-meson states are commonly parameterised in terms of 
the i?-parameters Bb and Bb^, for instance 

(^0|OLL|i?°) = ImlfBBB , (1.2) 

for £ = d. The operator Oll can be decomposed into parity-even and parity-odd compo- 
nents Ovv+AA and Ova+av- In the Standard Model only the parity-even part Oyv+AA 
makes a contribution to S-meson mixing. The S-parameters encode the long-distance 
effects of the strong interaction and must be determined in a non-perturbative approach 
such as lattice QCD. Indeed, various lattice estimates of Bb and Bb^ have been published 
by several authors in recent years [1~12]. 

It is well known that relativistic 6-quarks cannot be simulated directly for currently ac- 
cessible lattice spacings. Several formalisms for treating 6-quarks on the lattice, based on 
Heavy Quark Effective Theory (HQET) [13, 14], non-relativistic QCD [15], on-shell im- 
provement for relativistic quarks [16,17], as well as finite-size scaling techniques [18], have 
been developed and applied. Obviously, some, if not all, of these approaches imply certain 
approximations or assumptions whose validity and intrinsic accuracy must be investigated. 

In order to yield useful phenomenological information, operators like Oll must be renor- 
malised. If the regulator breaks chiral symmetry, as is the case for Wilson fermions, the 
renormalisation of OlL) which has a particular chiral structure, is complicated by the fact 
that it undergoes mixing with operators of different chiralities. Therefore, in addition to 
an overall logarithmically divergent, multiplicative renormalisation factor, one must also 
determine finite subtraction coefficients. 

The analogous case of K^ — K^ mixing, in which all fields that appear in the corresponding 
four-quark operator are treated relativistically, has been studied in ref. [19]. There the 
renormalisation and mixing patterns of a general set of four-quark operators were classified 
according to their transformation properties under certain symmetries. In particular, it 
was shown how the mixing due to explicit chiral symmetry breaking implied by the Wilson 
term could be isolated and absorbed into mixing coefficients. Another important result 
of [19] was the observation that the parity-odd component OvA+AV is protected against 



mixing by discrete symmetries. 

In this paper we adopt a similar strategy to extend the analysis of ref. [19] to the case 
where the 6-quark is treated at leading order in HQET, i.e. in the static approximation. In 
particular, we show how the heavy quark spin symmetry, in conjunction with transforma- 
tion properties under spatial rotations, as well as discrete symmetries like parity and time 
reversal can be used to constrain the renormalisation patterns of a general set of static-light 
four-quark operators. One key result is that it is possible to find a basis of parity-odd oper- 
ators that renormalise purely multiplicatively. This allows to devise a strategy aimed at a 
non-perturbative determination of the renormalisation factors required for the calculation 
of S-parameters, without the need to determine finite subtractions. 

To this end we use twisted mass QCD (tmQCD) [20] as a discretisation for the light quark 
fields, which allows us to map parity-even operators to parity-odd ones. By employing 
the Schrodinger functional (SF) [21], the anomalous dimension of the latter can then be 
determined non-perturbatively in complete analogy to the case of K^ — K^ mixing studied 
previously in [22]. Thus, in order to reconstruct the phenomenologically relevant matrix ele- 
ment of OlLi one only needs to determine the renormalisation properties of multiplicatively 
renormalisable operators, even in regularisations that break chiral symmetry explicitely. 

In addition to explaining how the renormalisation properties of static-light four-quark op- 
erators can be constrained, another purpose of this paper is to identify - in the spirit 
of [23] - suitable finite-volume renormalisation schemes based on the SF, to be used in a 
forthcoming non-perturbative calculation. To this end we have computed the anomalous 
dimension in perturbation theory at NLO for the complete basis of four-quark operators 
in several SF schemes. 

Of course, the complicated mixing patterns one is confronted with when using Wilson 
fermions can be avoided by using discretisations for the light quarks which obey the 
Ginsparg- Wilson relation. First steps in this direction have been taken in ref. [24, 25]. 
However, in this work we show, by using symmetry properties and tniQCD, that the renor- 
malisation of static-light four-quark operators describing B^ — B^ mixing can be studied 
in an equally simple framework for Wilson-type regularisations. Non-perturbative renor- 
malisation can thus be implemented in a straightforward manner and at much reduced 
computational cost. 

This paper is organised as follows: in section ^ we discuss how the transformation properties 
under various symmetries constrain the mixing patterns of static-light four-quark opera- 
tors. In section y we formulate a set of renormalisation conditions for the operator basis 
within the Schrodinger functional. Section ^ describes the perturbative calculation which 
yields the NLO anomalous dimensions of the operators for a set of Schrodinger functional 
renormalisation schemes. In section g we discuss the use of tmQCD to compute the physical 
matrix elements for B^-B^ mixing using multiplicatively renormalisable operators. Our 
conclusions are presented in section ^. Technical details regarding the use of symmetries to 
constrain the renormalisation pattern and the evaluation of lattice integrals are relegated 
to Appendices ^ and |^, respectively. Tables listing the finite parts of renormalisation 



constants and the NLO anomalous dimensions can be found in Appendix ^. 

2. Mixing of heavy-light four-quark operators in the static approximation 

In this section we study the mixing of AB = 2 heavy-Ught four-quark operators in which the 
heavy quarks are treated in the static approximation of HQET. Thus they are represented 
by a pair of static fields {iph^i^h), propagating forward and backward in time, respectively; 
their dynamics is governed by the Eichten-Hill action [26] 

S^'-'^'i^h^A] = a''Y,[Mx)V*oMx)-M^)VoM^)] , (2.1) 

X 

where the forward and backward covariant derivatives Vq, Vq are defined by 

VoV'h(3;) = - [Uo{x)i^h{x + aO) - i^hi^)] > 

V*oM^) = I [Mx) - Uoix - a6)-^Mx - aO)] . (2.2) 

The field iphi'iph) can be thought of as the annihilator( creator) of a heavy quark. Simi- 
larly, iphiiph) creates(annihilates) a heavy antiquark. Each field is represented by a four- 
component Dirac vector, yet only half of the components play a dynamical role, owing to 
the static projection constraints 

p+^Ph = V'h , i^hP+ = i^h , P+ = -(1 + 7o) ; 

P-A = ^-h, i^-hP- = h, P- = ^(l-7o). (2.3) 

Instead of the link variables that appear in eq. ( |2.2| ) one can consider more general defini- 
tions of the parallel transporter which enters the covariant derivative. A set of alternative 
discretisations was studied in [27], where it was found that adequate choices of parallel 
transporter lead to much improved signal-to-noise ratios in actual simulations. 

The light (relativistic) quarks are instead taken to be Wilson fermions, using either the 
plain Wilson action or its 0(a) improved version with a Sheikholeslami-Wohlert term [28]. 
The explicit chiral symmetry breaking induced by the Wilson term causes the mixing of 
operator of different naive chirality even in the chiral limit. 

We consider a complete basis of AB = 2 heavy-light four-quark operators which, for 
the sake of definiteness, we chose to contain two static fields "^/j and tpj^ while, in the light 
sector, we consider massless fermions with two distinct flavours -01 and ■02- We introduce 
a generic AB = 2 operator via 

^Sr^ = I [(V5hriVi)(v5;ir2^2) ± (V^hriV'2)(4r2V'i)] , (2.4) 

where ri^2 are Dirac matrices, and we adopt the notation 

^riT2 ± r3r4 = ^rir2 ^ ^r3r4 ■ (^•^) 



The complete basis of Lorentz invariant operators is given by the set of 16 operators 
parity-even: Qf = Oyv+AA ' parity-odd: Qf = O^a+av ' 

Q2 = ^SS+PP ' ^2 = C'gp+pg 

^3 = ^VV-AA ' ^3 = ^VA-AV ' 

Q4 = ^SS-PP ' ^4 — ^SP-PS ' (^•*^) 

which we have grouped according to their transformation properties under parity. Here 
V = 7^ , A = 7^75 , S* = 1 , P = 75 , and an imphcit summation over pairs of Lorentz indices 
is understood. We incidentally remind the reader that tensor structures like T = a^y or 
T = CTfip^i^ produce redundant operators in the static limit, due the projection constraints 

The description of IS.B = 2 transitions in terms of the static approximation of HQET 
implies that the operator Oll of eq. ( |1.1| ) is related in some particular way to the operators 
listed in eq. ( |2.6| ). Owing to the heavy-quark spin symmetry, one finds that Oll must be 
matched to a linear combination of Oyv+AA ^^^ ^ss+pp t^^l' ^'^'-^ thus those two operators 
are of particular interest to the study of B^ — B^ mixing. 

The operator basis in eq. ( |2.6| ) renormalises, in full generality, via a 16 x 16 matrix Z, 
the form of which can be constrained through symmetry arguments. A systematic method 
to carry out this analysis is given by the following prescription: 

(i) Construct the matrices $fc that implement, at the level of the operator basis, a 
maximal set of independent symmetry transformations that leave the action invariant. 

(ii) Impose the constraints 

Z = $fcZ$^ \ \/k. (2.7) 

The solution Z to this system of equations displays the constrained form of the 
renormalisation matrix. 

In most cases the constraint imposed by a given symmetry on the renormalisation matrix 
Z can be easily found out, while in a few cases (namely heavy quark spin symmetry 
and H{2>) spatial rotations) an explicit construction of the corresponding <l>fc matrices is 
required. We leave the explanation of this procedure to Appendix ^ and we present here 
the list of symmetries that have been used and their effect in constraining the matrix Z. 

Flavour exchange symmetry S. S exchanges the two relativistic flavours ipi and "02- 
Operators with superscript ± are eigenvectors of <1>5 with eigenvalues ±1 respectively. S 
thus prevents the mixing between the -|- and — sectors. This reduces the renormalisation 
matrix Z to a block-diagonal form, with two 8x8 blocks. 

Parity. Mixing among operators with opposite parity is excluded, and the renormalisation 
matrix Z is reduced to a block-diagonal form, where four 4x4 blocks describe the mixing 
of the parity-even and parity-odd operators among themselves. 



Chiral symmetry. It is used in the same way as in ref. [19]. In the chiral hmit, the 
continuum relativistic quark action is invariant under the finite axial transformation: 



ipk -^ ij^tpk ; Tpk -^ itpk75 ■ 



(2. 



Under this transformation we obtain: 
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(2.9) 



From this one sees that, were chirality respected by the regulator, Q^ would mix only 
with Q2 , and Q^ only with Q^ (and similarly in the parity-odd sector). This is not 
the case with a Wilson regularisation, for which the structure of chiral multiplets must 
be restored by combining operators with different naive chiralities [30]. The restoration 
of chiral properties is achieved by introducing the mixing matrices A^,il . Once the 
subtracted operators Q"^ = (1 + A'^)Q'^ and Q"^ = (1 + iI^)Q^ with the correct chiral 
properties have been constructed, they will mix like in the continuum with renormalisation 
matrices Z ,Z . We choose the matrices Z ,Z ,A ,Jl such that: 
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(2.10) 
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(2.11) 



This choice is convenient because it is easy to show (for instance by using Ward identities) 
that A , il and the product Z (^Z )^^ all depend only on the bare coupling gO) while 
Z and Z alone contain the continuum-like dependence on the renormalisation scale. 

Heavy quark spin symmetry and -ff(3) spatial rotations. Further constraints can 
be obtained from the heavy quark spin symmetry and cubic rotations. The procedure is 
slightly involved and we leave its description to Appendix ^ It applies identically to both 
parity-even and parity-odd sectors, and below we provide the expressions for the latter — 
results for the parity-even sector are obtained by simply replacing the symbols Q, Z^ JX 
with Q, Z, A. After imposing the constraints Z = $fcZ$^ one finds that it is possible to 
rotate (^) into a new basis 



{Q^tMt.Q^tMtf 



(Qf,Qf +4Qf,Qf + 2Qf,Qf 



2Qf)^, 



(2.12) 



in which the scale-dependent mixing is completely disentangled (even though some scale- 
independent mixing remains): 



(2.13) 
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Time reversal. Up to this point, the renormalisation mixing of the parity-even and 
parity-odd sectors has the same matrix structure. We now consider the effect of a time 
reversal transformation of both the static and relativistic quark fields. To that purpose, it 
is convenient to rewrite (|2.4D in the form 



Orir^ = 9 [(^i^+riVi) [^P-T2^2) ± (^^-riV'2) (^P+raVi)] 



with 



It is then easy to apply the time reversal transformation 



^(x) -^ *(x^)757o, 



1,2 



(2.14) 
(2.15) 

(2.16) 



where {xo,x.Y = (— xo,x). 

The transformation rules for the operators (apart from the r reflection of the space- 
time coordinates) are easily found to be 



Qf ^ ±Qf 



Qf ^ TQf 



Qt- 


-±Qf 


Qi- 


-.±Qf 


Qi- 


-±Qi 



Qf- 


-TQf 


Qf- 


-^±Qf 


Qf- 


-^±Qi 



(2.17) 



and identical ones for the {Q' , Q') basis. It is then clear that no new constraints arise for 
the parity-even operators and for Z'^ ^ while the invariance of the scale independent mixing 
under time reversal immediately yields 



R 



i± 



0. 



(2.18) 



This proves the purely multiplicative renormalisability of the operator basis ( 2.12| ). As a 
note of reference, we would like to point out that the absence of mixing is already manifest 
at the one-loop perturbative level in eqs. (20)-(23) of [6]: it is enough to project both sides 
onto parity eigenstates and change to the basis in eq. ( |2.12D to find that mixing in the 
parity-odd sector is absent at one loop. 



Finally, we can rotate back to the standard basis (| 
for the renormalisation matrices: 



I), obtaining the following form 
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(2.19) 



(2.20) 



il^ 



(2.21) 

(2.22) 



For later convenience, we denote by A the matrix responsible for the change of basis (|2.12|), 



i'± 



such that Qj = AjjQ. and Q^ = AijQ- . Equation (|2.19D reproduces the result of [24]. 



Now we will pursue a strategy to compute the renormalisation matrices non perturbati- 
vely using Schrodinger functional techniques. For parity-even operators, the determination 
of the subtraction coefficients in eq. ( p. 20 ) can be achieved by implementing suitable axial 
Ward identities. However, as we will show later on, the use of tmQCD techniques allows to 
obtain all the AB = 2 physical amplitudes of interest in the Standard Model from matrix 
elements of the operators Qf , Q^. Therefore, from now on we will concentrate exclusively 
on the renormalisation of Q^ . 



3. Renormalisation conditions in the Schrodinger functional 

Having constrained their renormalisation patterns by imposing the symmetries of the the- 
ory, we now proceed to specifying suitable renormalisation conditions on the four-quark 
operators ( 2.12| ). To this end we will consider a set of correlation functions defined in 
the Schrodinger functional (SF) which will serve to impose the renormalisation conditions 
at the non-perturbative level. The SF formalism [21], which was initially developed to 
produce a precise determination of the running coupling [31-33], has been extended to 
various other phenomenological contexts. These include the study of quark masses [34-36] 
and decay constants [37-39], the computation of moments of structure functions [40,41], 
the Kaon i?-parameter [22,23,42], and the static-light axial current [43,44]. The reader 
is referred to [45] for detailed explanations of the framework and the standard notation, 
which are not repeated here. 

Our construction of the relevant SF correlators extends the one described in [23] to 
the static case. We start by introducing the following bilinear boundary source operators. 



projected to zero external momentum, 



x,y 



Ks^r]=a^^CU^)TC,iy), (3.1) 



x,y 

where F is a Dirac matrix, the flavour indices si,2 can assume both relativistic and static 
values, and the fields ( and C' represent functional derivatives with respect to the fermionic 
boundary fields of the SF, see [45]. The choices for F are limited by the boundary conditions 
imposed on quark fields. At the level of boundary quark and antiquark fields they imply 

C(x) = P^a^) , C(x) = C(x)P+ , (3.2) 

with P± = 1(1 lb 7o), and similarly for C, (^'. Therefore, in order to have a non-vanishing 
source field Ss^sji^], F must anticommute with 70."^ 

Starting from the bilinears in ( p.l[ ), we have to construct suitable boundary sources to 
probe four-quark operators. In order to define a set of non-zero correlators in the massless 
theory, which can then be used to impose renormalisation conditions in the chiral limit, the 
probe should be parity-odd. We further require it to be invariant under the group H{3) of 
lattice rotations in three dimensions. This leads us to introduce a third "spectator" light 
quark and consider, as the simplest possible choice, a generalised boundary source made 
of three bilinears, 

W[Fi,F2,F3] = Si^[Fi]S23[F2]S'3;^[F3] , (3.3) 

two of which are localised at the boundary xq = 0, while the third lies at the other bound- 
ary, namely xq = T. Odd parity and rotational invariance are then assured through an 
appropriate choice of the Dirac structures [Fi, F2, F3], and a maximal set of corresponding 
probes is given by 

5«=W[75,75,75] , (3.4) 

1 ^ 

'^^'^ = 6 ^ ^kimm7k,7inm\ , (3.5) 

k,l,m=l 
1 ^ 

5(3) = _^>V[75,7,,7fc], (3.6) 

fc=i 

1 ^ 

S^'^ = i;Y.^hk,i5nk] , (3.7) 

-^ k=i 

1 ^ 

5(5) = ^>V[7fe,7,,75]. (3.8) 



3 

fc=i 



^Allowing for a non-vanishing angular niomentum would relax this constraint, but, since it would most 
likely lead also to worse signal-to-noise ratios in numerical simulations, we do not pursue this approach. 




Figure 1: Diagrammatic representation of correlation functions: the four-quark correlator 
•^fe (-^o) (left), the boundary-to-boundary static-light correlators fi.hi,ki^hi (center) and the 
boundary-to-boundary light-light correlators /i,zi,fci,iz (right). Euclidean time goes from left to 
right. The double blob indicates the four-quark operator insertion. Single lines represent relativis- 
tic quarks and double lines denote the static ones. 



The four-quark operators can now be treated as local insertions in the bulk of the SF, and 
their correlators are naturally defined as 



r, 



■±;(«), 



xq) 



L-3(Q'±(x)5(^)), . = 1,...,5. 



(3.9) 



A pictorial interpretation of (|3.9| ) is provided by the left diagram of Figure |l|. 

It has to be observed that, due to the symmetries of the static approximation for 
the heavy quarks, not all of the above correlation functions are independent. By using 
the explicit spin structure of the static propagator, some straightforward though tedious 
algebra leads to the constraints 
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(3.10) 



These relations show that only 24 out of the 40 correlation functions defined in eq. ( p. 9] ) 
are independent. We stress that this result is exact; as a cross-check, later on we will find 
these identities to be explicitly verified at one loop in perturbation theory. 

In order to isolate from eq. ( |3.9| ) the ultraviolet divergences that are due to the bulk 
operator and absorb them into a renormalisation factor, one has to address the renormal- 
isation of the boundary fields. The ultraviolet divergences of the latter are cancelled by 
defining suitable ratios of correlators for which the renormalisation factors of the boundary 
fields drop out. To this end we introduce a set of boundary-to-boundary light-light and 
static-light correlators, 

^ -\(D\^{l^\0^x{l^\) , (3.11) 



fi 



hi 



fll 
/l 



2L6 

1 
'2I6 



(0'l2[75]02l[75] 



(3.12) 



k=l 



(3.13) 



whose valence structure is represented as well in the middle and right diagrams of Figure |l[ ^ 
A suitable combination of such correlators must comprise the same number of static and 
light boundary fields as (|3.9| ) . The simplest examples are given by the ratios 






(3.14) 



involving either the correlator /j' or k'^. However, it is easy to generalise these conditions 
by introducing an arbitrary real parameter a. Hence, we consider the following ratios of 
correlation functions: 



fhyUjl/2-a^j^ll^ 



K;ai^o)= ..ulu^Zu.. - (3-15) 



Although there is no real a priori restriction on the value of a, it is clear that "natural" 
values should lie in the interval [0,^]. This freedom, together with the choice of the 
boundary source and the 6'-angle of the SF, can be used in a later stage to tune the 
optimal renormalisation schemes, with the aim of having small NLO coefficients in the 
corresponding anomalous dimensions. This is important in order to control the systematics 
of the perturbative matching to continuum schemes at high energy scales. 

For the moment we observe that the ratios (p. 15) are free of boundary divergences. 



and consequently we impose the renormalisation condition, 

2S^'^hi;!^(T/2) = hi;!^^(T m\,o=o , (3.16) 

where all the correlation functions are computed in the chiral limit. This fixes non- 
perturbatively the renormalisation constant 2^,.^ at the scale // = 1/L. As usual, the 
Z' factors depend upon every calculational detail with the only exception of the leading 
log, which is universal. In order to operatively define a renormalisation scheme, a com- 
plete specification of the parameters that concur to quantify ( 3.15| ) is required. We briefly 
summarise them: 

• the possible presence of a Sheikholeslami-Wohlert (SW) term in the lattice action for 
the light quarks; 



the choice of the gauge parallel transporter in the covariant derivatives, eq. (p.2[), in 
the static action; 

the Dirac structure of the boundary source; 

the value of the angle 6 entering the spatial boundary conditions of the SF; 

the value of the parameter a in eq. ( |3.15| ); 

the ratio T/L between the time and the spatial extension of the SF. 



^The static-light counterpart of fc" is not considered as, due to the symmetries of the static limit, it is 
identical to fiK 
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The last four conditions fix the renormahsation scheme, while the first two only introduce a 
regularisation dependence in the renormalisation constants. We stress at this point that the 
running of the operators is a continuum property, i.e. it is independent of the discretisation 
chosen. The latter only affects the way in which the continuum limit is approached, and - 
in the case of the action for static quark fields - the signal-to-noise ratio in the simulations. 
In this paper we consider the Wilson action (with and without SW term) for the light 
quarks, and the Eichten-Hill action for the static ones. As for the parameters that fix 
the renormalisation scheme, we will consider three values of (namely 9 = 0.0, 0.5, 1.0) 
and two values of a (namely a = 0.0, 0.5). Furthermore, we fix T = L. Taken together 
with the possible independent choices of boundary sources, this leaves us with 12 different 
renormalisation schemes for the operators Q\ and Qg , and 24 schemes for Qg ^^'^ Q'4 ■ 
The number of independent renormalisation conditions is twice this figure, as we consider 
two different actions. 

4. A perturbative study 

We now proceed to studying the renormalisation of the operators Qi,...,QJ at one- 
loop in perturbation theory. These are the operators that will enter AB = 2 effective 
Hamiltonians in the static limit. The purpose of this study is threefold. First, it provides 
an explicit check of the expected mixing pattern. Second, it will allow us to compute 
the NLO anomalous dimension of the operators in the SF schemes defined above; this is 
done via a standard one-loop matching procedure to continuum schemes where the NLO 
anomalous dimension is already known. Finally, we can work out lattice artefacts of the 
step scaling functions in one-loop perturbation theory. They can then be compared to 
the ones in the relativistic case discussed in [23], in order to obtain information about the 
expected size of discretisation effects for quantities involving static fields. In principle, the 
discretisation effects computed in perturbation theory could subsequently be subtracted 
by hand from the non-perturbative Monte Carlo data, with the aim of exerting a better 
control of their continuum extrapolation, as pursued in [22,46,47]. 

4.1 Scheme dependence of NLO anomalous dimensions 

In order to proceed, some notation has to be fixed. The scale dependence of operator 
insertions in renormalised correlation functions is described by the RG equation 



^'^ (^|l + ^1 + ^'^§x + ™^) - ^{ 



(Q|)r = , (4.1) 



-J 



where /3 is the /3-function for the coupling, r is the anomalous dimension of the quark 
mass, and 7"*" is the operator anomalous dimension matrix. We have also included a term 
which takes into account the dependence on the gauge parameter A in covariant gauges, 
characterised by the RG function f3\, defined as 

.^ = AA. (4.2) 
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This term is absent in schemes hke MS (irrespective of the regularisation prescription) or 
the SF schemes introduced in section ^, but is present e.g. in regularisation-independent 
(RI) schemes, which wih be considered later on. If we choose to work with the operator 
basis in eq. ( 2.121 ), the matrix structure of the RG-equation simpUfies, and the evolution 



of the various operators is determined by a set of scalar anomalous dimensions. 

In what follows, we focus on mass-independent renormalisation schemes, for which 
the RG functions depend only upon the coupling. We take the following form for their 
perturbative expansions: 



oo oo 

J2k 



Pig) = -/ E ^f^a'' , Pxig) = -/ E ^'^g' 

fc=l A;=0 

oo oo 

r{g) =-g'Y.dk9'' , itM = -g'T.^^9'' , (4.3) 



fc=0 fc=0 



with universal coefficients 



Oo 



1 fll 2 1 , I { N 

N--N,\ , do = __J3 



(47r)2 [3 3 'J ' ^ (47r)2 [ N 



t2 






(47r)4 [3 V 3 

The LO coefficient of the anomalous dimension matrix j~^'^^> is universal as well, and 
has been calculated in [29] for the first operator of the basis and [6] for the rest. The 
non- vanishing elements of the anomalous dimension matrix read 

-;{0) 1 /^.. 3 






(47r)2 V N 



^ir = -7£w ■ ("-9) 



-;{0) _ 3/2 
■ (47r)2 



74"/"^ = -TTT^I ' (4-10) 



-;(0) 1 /„,. 6 



^" =-(4^r-^i- <"'' 



A covariant rotation of this matrix to the diagonal operator basis ( 2.12| ), i.e. 7'^' 



;(o) 



A.j'^'^WA ^, gives the LO coefficients of the multiplicatively renormalisable operators, 
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namely 



'+;(o) 

7i 



(47r)2 



3N 



3_ 

N 



/+;(0) 
72 



'+;(0) 
73 



3iV + 3-A 



'+;(o) 

74 



(47r)^ 



1 



37V -4 



7_ 

N 



3N 



"^ iV 



(4^)2 \^"-' ' " Nj ' '4 (47r)2 V"^' " Nj ■ ^^"-^^^ 

By contrast, the NLO coefficient is scheme-dependent. The perturbative matching 
procedure that ahows to express its value in the SF scheme in terms of the value in a 
reference scheme has been derived in [48] for the case of multiplicatively renormalisable 
operators. The formalism can be trivially extended to situations where mixing occurs and 
a gauge covariant reference scheme is assumed. The renormalised operators and coupling 
constant are first related in the two schemes through a finite renormalisation, 

9SF = '^(;(5rcf)5'ref ' 
(Q+Sp)r = <^+(5ref)(Qicf)R • (4-13) 

The matching coefficients X are then expanded in powers of the coupling constant, 



k=l 



(4.14) 



and the requirement of formal invariance of the RG-equation under a change of renormal- 
isation scheme leads to the two-loop matching relation 

^3+f ) = ^+f ) + [^+;(i),^+;W] + 26oA'+^(i) + b'oX§^X^'^'^ - I'-^^'^X!^^) , (4.15) 

where the symbol [•,•], which is absent in the case of multiplicative renormalisation, repre- 
sents the ordinary matrix commutator. It should be stressed that the choice of the reference 
scheme is irrelevant. In fact, a good consistency check on the result for 7gp is provided 



by computing the RHS of ( 4.15| ) for several different reference schemes. 



Finally, we point out that the lattice is currently the only known regularisation of the 
SF, for which perturbative calculations of fermionic observables can be operatively per- 
formed.^ If the reference scheme is defined in the continuum, the operator matching must 
take into account both a change of regularisation and a change of subtraction prescription. 
Accordingly, Xq ' must be computed as the difference of two matching coefficients to an 
intermediate scheme, namely 



A,^ 



;(i) 



;(i) 



;(i) 



'•SF.ref ~ '^SF.lat '^ref,lat • (4.16) 

The "lat" scheme is by definition the minimal subtraction lattice scheme, where the renor- 
malisation constants are polynomials in ln(a/x) without finite parts. Consequently, X^^\^^, 
which provides the matching between SF and "lat" , can be obtained from a one-loop calcu- 
lation of the renormalisation constant in the SF scheme with a lattice regularisation. The 
matching coefficient -^j-ef lat between the reference scheme and the lattice can be instead 
retrieved from the literature for som e choice of the reference scheme, such as MS or RI. 

^A recent proposal to perform the matching directly in dimensional regularisation has been presented 
in [49]. 
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4.2 Perturbative expansion of SF correlation functions 

We now describe the one-loop calculation of the SF renormalisation constants introduced 
in section y. The perturbative procedure is fairly conventional, and we include it just for 
completeness. We start by expanding all the correlation functions previously introduced 
in powers of the bare coupling, 



X = X<»' + gl 



orriQ 



c 



+ 0(gt), (4.17) 



where X is one of J-j^ , /{* , f[, k{, or a linear combination thereof. The derivative term 
in square brackets is required in order to set the correlation function X to zero renormalised 
quark mass, when each contribution to the RHS is calculated at zero bare quark mass, as it 
will be assumed. As for the numerical value of rric , we use the numbers provided by [50], 
i.e. 

.,. f -0.20255651209 Cf (csw = 1), ^ iV" - 1 ,,,,, 

amy-' = { ' ' Cf = r^, (4.18) 

[ -0.32571411742 Cf (csw = 0), 2A^ 

The SF renormalisation constants, defined in (|3.16| ), admit an analogous expansion, 

4+f )(50,a/L) = l + glz'^-;^^'^'\L/a)+0{gt) . (4.19) 

The explicit expression of the one-loop order coefficient Z^,.^ in terms of the perturbative 
expansion of the four-quark and the boundary-to-boundary correlators can be obtained by 
inserting ( 4.17 ) and (|4.19| ) into the renormalisation condition (|3.16| ). One then obtains 






y-^k -^k -^k ) xo=T/2 

Ml) f«(i) ^(i);^f"(o)1 




Contributions containing the subscript "6" arise from the boundary terms that are required 
in addition to the SW term in order to achieve full 0(a)-improvement of the action in the 
SF [45] . Obviously, these contributions are not present when the unimproved Wilson action 
is chosen for the light quarks. From now on we will set them to zero also when the action 
is 0(a) improved, as they will not affect the continuum limit extrapolations involved in 
the computation of NLO anomalous dimension, and their contribution to cutoff effects is 
negligible. 
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The evaluation of the RHS of ( 4.20 ) requires the calculation of the Feynnian diagrams 
depicted in Figures ^ and ^. The one-loop expansion of the boundary-to-boundary cor- 
relators fi and k^i is known from [51], while /^' has been studied perturbatively in [43]. 
Accordingly, the only new diagrams which need to be calculated are the ones that contribute 
to the one-loop order coefficient of the four-quark correlators. Two groups of diagrams can 
be identified: the self-energies correct the valence fermion propagators through a gluon 
emission with subsequent absorption by the same leg, and the vertex diagrams correct the 
operator insertions through the exchange of a gluon between two legs. Each of them can 
be expressed as a loop sum of a Dirac trace in time-momentum representation, where the 
spatial coordinates are Fourier transformed. These sums have been performed numerically 
in double precision arithmetics using a C++ code, for all the even lattice sizes ranging 
from L/a = 6 to L/a = 48. The results have been checked by an independent Fortran 
90 program, also in double precision arithmetics. The behaviour of the renormalisation 
constants thus obtained, as functions of the lattice size L/a, is expected to conform to the 
standard asymptotic expansion 

oo 

^^^'''^iLH = E ilYH^ + <ML/a)} , (4.21) 

which can be used in order to extract the universal LO anomalous dimensions and the finite 
constants peculiar to the schemes, that is to say, the coefficients Sq and r^, respectively. 
The latter represents the matching coefficient introduced in ( |4.16 ) in the diagonal basis, 
namely 

<SF,iat = <o • (4-22) 

An efficient numerical technique to isolate these coefficients, based on a blocking procedure 
of the function at neighbour lattice sizes, has been introduced in [52]. Details about its 
application to the case at hand are provided to Appendix |^. Numerical values of the 
coefficients r^^ for the various schemes introduced in section |3| are reported in Tables ^ - 

& 

4.3 Matching to continuum schemes and consistency checks 



The NLO anomalous dimension matrix of the operators ( p.6D in continuum schemes can 
be found in [6], together with the one-loop matching relations to the minimal subtraction 
lattice scheme. The regularisations employed in [6] are DRED and NDR, and two possible 
subtraction prescriptions are considered, namely MS and RI. An attractive feature of the 
latter is the independence of the corresponding anomalous dimension from the choice of 
evanescent operators (EO), which complicate the mixing pattern in d = 4 — 2e dimensions. 
As a consequence, it is trivial to perform a rotation of the anomalous dimension matrix 
in the RI scheme to a different basis of the physical operators, such as ( |2.12| ), without the 
need to address subtleties related to the definition of evanescent contributions. The choice 
of RI as a reference scheme is therefore convenient in order to make use of the two-loop 
matching relation ( 4.15| ) in the diagonal basis ( p. 12 ). 
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Figure 2: Feynman diagrams of the self-energy type. 








Figure 3: Feynman diagrams of the vertex type. 
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Results reported in [6] refer to a perturbative expansion in powers of the MS-coupling. 
We therefore need the matching coefficient in eq. ( 4.13|) , which relates gsF to fl^, to 
one-loop order. This has been calculated in [53] and is given by 



X^^^ = 26o ln(/iL) - i- (ci,o + ci,iiVf) , 
ci,o = 1.25563(4) , ci,i = 0.039863(2) . 



(4.23) 



The NLO anomalous dimension of the operator basis ( |2.12| ) in the Feynman gauge (A = 1) 
and NDR regularisation^ , obtained from the covariant rotation 7'^^-* = A7'^)A~^, is a 
diagonal matrix whose non-zero coefficients read 



7l;RI 



(4^) 



131-F87r2 2 „,. l-27rV 30 -h 47r^ 57 + 167r^ 



12 



-N^ + 6N 



+ 



3iV 



12Af2 



, 5 5 

+ iVf -N 

'3 3N 



72;RI 



(4^) 



73;RI 



(4vr) 



74;RI 



(4^) 



131 + Svr^ 



12 



9 

,,,.,, 40 55 ^ 

+ Nf{-N 

* ' 3 9 9A^ 



301 + 67r2 18 - in^ 87 - Wn^ 



9 



+ 



3N 



+ 



12A^2 



131 + 87r2 2 83 ,_ 309 + 87r^ 24 - 47r^ , 21 - 47r^ 



12 



-N^ N + 

6 



12 



3A^ 



+ 



3iV2 



,5 10 5 

+ Nf{-N + 

'33 N 



131 + 8tt' 



12 



N^ + —N + 
6 

'3 3 N 



309 + 8tt^ 42 - 47r2 33 - Att^ 



12 



+ 



3A^ 



+ 



3iV2 



(4.24) 



The same rotation can be applied to the one-loop operator matching matrix X^-^ ■^^^. In this 
case the analytic dependence upon the gauge parameter A is needed in order to account 
for the derivative term included in the two-loop matching relation ( [4.15| ). With A^ = 3, 



^Although the four-quark operators are renormalised according to the RI scheme, which is independent 
of the regularisation prescription, the strong interaction Lagrangian is renormalised in MS . This introduces 
a spurious dependence of the NLO anomalous dimension upon the choice of the regulator. 
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one has 



X\ 



'(1) 

l;RI,lat(wilson) 



x: 



'(1) 



2;RI,lat{wilson) 



x: 



7(1) 



x\ 



3;RI,lat(wilson) 
'(1) 



4;RI,lat(wilson) 



1 

(47r)2 


1 

(47r)2 


1 

(47r)2 


1 

(47r)2 



10 8 



A + [Dll - Drr) 



3 3 ^^ 



10 

T 



;A + - {2Dlr + 2DfK - 45Ri - Df,L) 



Y~l^^\ (2^i^ + 2^LH + 45i?L + ^II) 



(4.25) 



for light quarks regularised with the pure Wilson action. If the 0(a) improved action is 
used instead, one has to add to them the matching factors between the two actions, viz. 



X, 



X, 



'(1) 

1 ; lat (sw) ,lat ( wilson) 

'(1) 

2;lat(sw), lat (wilson) 



1 



X. 



'(1) 



3;lat(sw), lat (wilson) 



X, 



'(1) 



4;lat(sw), lat (wilson) 



(4^)2 
1 

(4^ 

1 
(4^ 

1 
(4^ 



4,7 1 / 4 , 
3'' 3 3 



-V-V 

3'' 9 



0.038033(2) 
0.040240(2) 
0.034253(2) 
0.037720(2) 



(4.26) 



where the lattice integrals /^, v^ and w^ are discussed in Appendix y. Numerical values of 
the D-coefficients, expressed in [6] as linear combinations of a basic set of lattice integrals, 
are reported in Table y of Appendix y, where a new computational method to improve 
their numerical accuracy is also described. The factors in eq. (|4.26| ) are obtained from the 
coefficients denoted D in [6], after subtracting the contributions coming from the 0(a) 
improvement of the four-fermion operators. 



All the ingredients needed to evaluate the RHS of eq. ( 4.151) have now been specified. 



The absence of operator mixing in the diagonal basis ( 2.12 ) implies that the commutator 
term in eq. ( 4.15| ) is identically zero. NLO anomalous dimensions in the previously intro- 
duced SF schemes follow from a straightforward use of eqs. ( [4. 121) , ( [4.16 ) and ( 4.22| )-( 4.25| ). 



We have collected the ratios of T^.gp to the corresponding LO coefficients 7^ in Tables 



10-15. In the matching we have employed the values of rtn obtained with the pure Wilson 



action for light quarks, as they tend to display a better behaved continuum extrapolation, 
after the 0(a) contributions have been removed through blocking. 

In order to check our results, we have also derived the SF NLO anomalous dimensions 
using MS as a reference scheme. The matching procedure, rather delicate in this case, 
must take into account the role played by the EO in fixing the finite contributions to 
the NLO anomalous dimension matrix 7^ . A naive rotation of the latter is potentially 
hazardous without reconsidering the choice of the EO. An alternative approach is to work 
within the original basis (p.6|), to which the results in [6] refer, and then rotate the one-loop 
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matching coefficients from the SF to the lattice scheme according to the inverse rotation 
'^SF lat ~ ^"^"^SF lat^- This is certainly possible, as the computation of such coefficients 
is performed on the lattice in d = 4 dimensions, where no EO contributes. Of course, 
the commutator term in eq. ( [4.15 ) must be included in this case, whilst the gauge term 



proportional to 6q is not present. Once the NLO anomalous dimension matrix has been 
obtained in the SF scheme, a straight rotation back to the diagonal basis (2.12|) yields the 



,/+;(!) 



scalar coefficients T/j.gj? • This procedure has been applied using either DRED or NDR 
regularisations. In both cases we obtain the same results as with RI in the diagonal basis. 

We have also verified that the difference between the finite parts of the SF renormal- 
isation constants with improved and unimproved Wilson light quarks coincides with the 
values in eq. ( [4.26| ). The numerical values of these finite matching constants are indeed in 
perfect agreement with the analogous SF quantities. 

We finally concentrate on the numerical values of the NLO anomalous dimension co- 
efficients in the SF. A comparison with the case studied in [23], where the four-fermion 
operators contain only relativistic quark fields, shows that in the present case the variation 
of the anomalous dimension due to different choices of the SF boundary sources in the 
renormalisation condition is much less pronounced. Also, the non-perturbative identities 
in eq. ( p. 10 ) are verified explicitly by the one-loop results. The dependence on the value of 



the parameter a is very small, too. Finally, at 6 = 0.5, which is commonly employed in non- 
perturbative studies of SF renormalisation, the values obtained for the ratios T^-gj? /T^ 
are relatively small, pointing towards a good convergence of the perturbative series, save 
for Q2 7 where they are close to —0.5. The question whether this is a relevant source of 
uncertainty in the NLO matching of renormalised matrix elements to continuum schemes 
is left for future studies. 

4.4 One-loop order cutoff effects in step-scaling functions 



The non-perturbative RG-evolution of the four-quark operators in the diagonal basis ( 2.12| ) 
is obtained through the computation of the step-scaling functions 






un. ai ^Lj I 

(4.27) 



g^{L)=u 



These ratios of renormalisation constants provide the operator running between the scales 
[i = 1/L and fi = 1/2L. The advantage of introducing such ratios is related to the com- 
pensation of logarithmic divergences between numerator and denominator, thus resulting 
in a finite continuum limit. Cutoff effects can be therefore completely decoupled from the 
continuum RG-evolution. We are concerned here with the perturbative expansion of ( 4.27] ) 
in the renormalised coupling, that is 

+;{«)/ N 1 I +;(s.i) i +;(s,2) 2 , r\^ 3\ ia r,o^ 

The first two terms of this expansion depend upon the LO and NLO anomalous dimension 
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1,0 running 

NT,0 running at 0- 0.0, a- 0.0 

NLO running at 6= 0.5, a= 0.0 

NLO running at 9= 1.0, a= 0.0 




[,0 runn 

NLO runn 
NLO runn 
NLO runni 



ng at 8- 0.0, a- 0.0 
at 6= 0.5, a= 0,0 
ng at 8= 1-0, a= Q.O 




Figure 4: On the left (right) side the step scahng function of Qi^(Q2^) at NLO and A^f = is 
plotted vs. the squared renornialised coupling in the SF scheme. The boundary sources choice is 
s — 1^ and the a-parameter is set to zero. 



coefficients. They read explicitely 



a 



+;(^,i) 



k\a 



a 



;(«,2) 



k;a 



/+;{0) 



In 2 , 



In 2 + 



1/ '+;{0)n2 , u '+;(o) 



(In 2)2. 



(4.29) 



A graphical representation of ( 4.28| ) for the whole operator basis in some particular cases is 
provided by the four plots of Figures |^ and ^, in the range of values of ggp used in previous 
non-perturbative studies by the ALPHA Collaboration. 

The rate of convergence of the step-scaling functions toward the continuum limit at LO 
can be expressed in terms of the first non-trivial coefficient of the perturbative expansion 
(analogous to ( 4.28 )) of Ti^!^ {u,a/L) via the ratio 



;W 



{a/L) 



.+;{s,i) 

'k\a 



{a/L) 



a 



+;(s,i) 

k;a 



a 



+;(s,i) 

k;a 



(4.30) 



'+;(s,i). 



'+;is,i). 



where 

^t;!:''\^/L) = Z;:;:'^"(2LM - Z'::r'\L/a) . (4.31) 

In order to compare the perturbative lattice artefacts ( 4.301 ) with the ones obtained from 
the corresponding non-perturbative Monte Carlo simulations, the same definition of the 
critical mass, based on the PCAC Ward identity, should be adopted. This point has been 
extensively explained in [23], where the numerical values of arric (L/a) from L/a = 6 to 
L/a = 32 have been provided (c/. Table 3 in that work). That discussion will not be 
repeated here. Since our codes ran up to L/a = 48, we are in the position to extend the 
aforementioned table to include the additional points. The new numbers are reported in 
Table |l[ 
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1,0 running 

NT,0 running at 9- 0.0, a- 0.0 

NLO running at 6= 0.5, a= 0.0 

NLO running at 9= 1.0, a= 0,0 




[,0 runn 

NLO runn 
NLO runn 
NLO runn 



ng at 8- 0.0, a- 0,0 
at 6= 0.5, a= 0,0 
ng at 8= 1-0, a= 0,0 




Figure 5: On the lcft(right) side the step scahng function of Q'^{Q'^) at NLO and A^f = is 
plotted vs. the renornialised couphng in the SF scheme. The boundary sources choice is s = 1, and 
the a-parameter is set to zero. 



L/a 


am^J''{L/a)\c,^=i/CF 


am^^\L/a)\c,^=Q/CF 


34 


-0.20255637783 


-0.32544080501 


36 


-0.20255639414 


-0.32547023220 


38 


-0.20255640819 


-0.32549515390 


40 


-0.20255642028 


-0.32551644442 


42 


-0.20255643068 


-0.32553477599 


44 


-0.20255643965 


-0.32555067224 


46 


-0.20255644740 


-0.32556454600 


48 


-0.20255645412 


-0.32557672619 



Table 1: The one-loop coefficients of the critical mass as obtained from the PC AC Ward Identity 
at finite lattice size. For the parameter choices made here, the convergence to the values at infinite 
lattice size is quadratic/cubic in {a/ L), for standard/0 (a) improved Wilson quarks. 



In practice, non-perturbative simulations based on the Eichten-Hill discretisation of 
the heavy quark fields should better be avoided, given the bad intrinsic signal-to-noise 
ratio (2.1, |2.2| ) [27]. Nevertheless, it is instructive to compute lattice artefacts in pertur- 
bation theory for the Eichten-Hill action, if only to check whether the use of static fields 
enhances lattice artefacts with respect to the purely relativistic case. A comparison be- 
tween static-light and light-light four-quark operators in a typical situation is shown in 
Figure^ All data refer to Qf, employing a renormalisation scheme in which the boundary 
sources have a Dirac structure [75,75,75] and where the normalisation of the four-quark 
correlator involves only f[K Relativistic data are taken from [23]. Taken at face value. 
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the plot leads to the conclusion that the light quark action is the main responsible for the 
presence of relatively large lattice artefacts: once this has been chosen, static-light and 
light-light four-quark operators come to be affected by cutoff effects of a similar size. 
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Figure 6: Comparison between cutoff effects of the step scaling function of the static-light and 
light- light versions of the four-quark operator Q'j^ . Light quarks are described in terms of the 
unimproved (W) or improved (SW) Wilson action. Static quarks are always described by the 
Eichtcn-Hill (EH) action. The critical mass has been obtained from the PCAC Ward Identity. 

5. tmQCD for B°-B^ mixing 

As stated above, our interest in the renormalisation of Q^ stems from the fact that the 
physical AB = 2 matrix elements involving the operators Qf,Q2i with ipi = 'tp2 = i^i 
(where in practice the flavour label l denotes either a d or an s quark), can be mapped 
onto matrix elements of Q^ , Q.2 computed in some suitable tmQCD regularisation. This 
reduces to a minimum the uncertainties related to operator renormalisation, which takes 
place essentially in the same way as if exact chiral symmetry were present. Now we will 
construct specific tmQCD regularisations which realize this mapping. This technique is a 
generalisation of the ones already developed for Bk [20,22]. 

We will work in the so-called "twisted basis" , and concentrate first on the case of the 
Bb parameter where ipi = d. Let us consider a quark doublet ^^ = {u,d), for which we 
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specify the tmQCD action 

5*-QCD[V,V^]=a^5^ {i;{x)[D„ + m + tfiTS5]i'{x)} , (5.1) 

where D^ is the usual Wilson-Dirac operator (with or without a SW term). The choice of 
action for the other relativistic quark flavours is immaterial to the argument. The equiv- 
alence of tmQCD to standard QCD has been first established in [20]. Given a multi-local 
gauge-invariant operator 0{xi, . . . ,x„), the equivalence amounts to the identity between 
renormalised correlation functions 

(Or(xi,... ,x„))(Mh,^o) = (CR(a^i>---,a;n))(mH,,MR) , (5.2) 

which holds in the regularised theory up to cutoff effects, and is exact in the continuum 
limit. In the above expression, the relation between the operators O and O is provided by 
the axial rotation of the quark fields which relates QCD to tmQCD, viz. 

ip = exp (ia'j^T^ /2) ip 
ip = ip exp [iaj^T /2j , 

where the twist angle a is defined in terms of the renormalised mass parameters of the 
tmQCD action as 

tan(a) = -^ , (5.4) 

and the physical renormalised quark mass Af^ is given by 



Mr = ^m| + nl . (5.5) 

Let us now consider the static-light four-quark operators Cyv+AA ^^'^ ^VA+AV with light 
flavours tpi = '4'2 = d. We observe that the rotation (|5.3|) implies 



^VV+AA = cos(a)0^v_^AA - isin(a)O^A+AV • (5-6) 

In particular, at a = 7r/2, which is known as the maximally twisted case, ( ^^ ) simplifies to 

^VV+AA = ~^^VA+AV • i^-^) 

Exactly the same property holds for the operator SS -|- PP, for which one finds again 

^SS+PP ~ ~^^SP+PS • (^•^) 

This demonstrates explicitly that the matrix elements of Cyv+AA ^^^ ^ss+pP' responsi- 
ble for the particle mixing in the SM and within the static approximation of QCD, can be 
obtained from a computation of the matrix element of OyA , ^y and Ogp i ps iii tmQCD at 
a = 7r/2. Since in mass independent schemes, such as the SF, all dependence of renormal- 
isation factors on the mass parameters drops out, it is clear that tmQCD does not spoil 
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the renormalisation pattern of the operator basis ( ^.12 ). In particular, the combinations 



Q!^ = Oya_,_^y and Q2 = ^ya+pc^ ~^ '^^sp+ps renormahse purely multiplicatively. 

In case one is interested in the B^-B^ mixing amplitude, it is enough to maximally 
twist a quark doublet which contains the s quark, e.g. (c, s), the action for which would 
read exactly as eq. (^]^), save for the eventual introduction of non-degenerate masses for 
the two quarks of the doublet, along the lines of [54]. The action for a twisted (c, s) doublet 
would then read 

^tmQCD^^^^] ^^4^ {V;(x)[Z?„ + m + i^75]V(x)} (5.9) 

X 

with m = diag(mc,TOs) and /^ = diag(//c) A's); and the constraint 

tan(Qj = = — ■ — . (5.10) 

"^s,R "T-c,R 



A potential shortcoming of eq. (|5.9| ) comes about in case it is taken as the action for dynam- 
ical c, s quarks, since M^^r 7^ Mc;k would then induce a phase in the fermion determinant. 
One may then consider more sophisticated chiral rotations that keep the determinant real, 
as in [55]. If c,s are kept quenched, or interpreted as valence quarks, no such subtlety 
arises. 

We conclude that the use of suitable tmQCD regularisations avoids the need of de- 
termining mixing coefficients for the renormalisation of the matrix elements entering the 
B^-B^ amplitude in the static approximation. For an alternative analysis of operator 
mixing using a different tmQCD regularisation, we refer the reader to ref. [56]. 

An additional advantage brought in by the use of maximally tmQCD is the automatic 
0(a) improvement of bare matrix elements of the above four-fermion operators. This 
property does not hold, on the other hand, for the renormalisation constants computed 
within the SF schemes discussed in previous sections. In order to obtain 0(o) improved 
renormalization constants one should use modified SF schemes, as proposed in refs. [57,58]. 

In order to show that the bare matrix elements are automatically 0(a) improved we 
extend the argument in Appendix A of ref. [59]. The first observation is that, at maximal 
twist, the only 0(a) counterterms to the static action are proportional to the dimension 
five operators 

tr(/x2) (V^ft^;, + Vi^iVh) , (tr(/x))2 {i>hi'h + i^hi^h) ■ (5-11) 

(In the simple case /x = ^r^ there is one single counterterm proportional to //^(V'hV'fe + 
i^h'^^h)-) These counterterms merely generate a shift of the static quark self-energy ^ . The 
second observation is that all the possible 0(a) (dimension seven) counterterms will have 
the same static field content of the original (dimension six) four-fermion operator, and will 
differ from it only by the addition of mass factors or derivatives. It is then possible to 
extend the symmetry V x Vd x (/i — > — ^) of the relativistic action [59] (where V is the 



^Moreover, they are obviously absent in the quenched approximation. 
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physical parity and Dd is defined in ref. [60]) to include static quarks. V and D^^ will now 
be defined to be 





Uo{x) - 


- ^o(x-) 




Uk{x) - 


-^ Ulix'' -ak) 




V'(x) - 


-^ iiQibni^ix'') 


V ■ < 


^(x) - 


-^ i){x'')T^-i^-iQi 


/ • ^ 


i^h{x) - 


-^ Mx"") 




i^h{x) - 


-^ Mx"") 




^hix) - 


- -Mxn 




. i^hix) - 


- -Mx"") 


{xo,^y 


= {xq,- 


-x). Using V X 



V, 



U^.{x) - 


-^ Ul{-x - afi) 


V'(x) - 


-^ —iip{—x) 


^(x) - 


-> —itp{—x) 


i^hix) - 


-^ i'hi-x) 


i^hix) - 


-^ i'hi-x) 


Aix) - 


-^ iphi-x) 


Mx) - 


-^ %i)h{-x) 



(5.12) 



T^d ^ il^ ^' ~f^) one immediately concludes that all 
the relevant dimension seven operators have opposite parity with respect to the dimension 
six ones. Using the same arguments of ref. [59], one then concludes that no 0(a) appears 
in the Symanzik expansion of the relevant correlation functions. 

6. Conclusions 

In this paper we have shown that the renormalisation problem of heavy-light four-quark 
operators in the static approximation can be tackled for Wilson-like fermions without the 
need to perform finite subtractions. 

Owing to the presence of static quark fields, the flavour switching symmetries, which 
in the relativistic case have proved so useful [19,61] for imposing constraints on the mixing, 
are very much reduced. However, this lack is compensated by the heavy quark symmetry, 
spatial rotations and a set of discrete symmetries, such as time reversal. The emerging 
renormalisation and mixing pattern is then quite similar to the relativistic theory: while 
chiral symmetry breaking generated by the Wilson term induces mixing among different 
chiralities of parity-even operators in the lattice regularised theory, such mixings are com- 
pletely absent in the parity-odd sector. 

Twisted-mass QCD can be used to relate the operator bases in the parity-even and 
parity-odd sectors also in the static approximation. In particular, we have shown how 
to do this for the operators that contribute to B^-B^ mixing in the Standard Model, 
using a maximal twist setup that brings in, as a bonus, the potential for automatic 0(a) 
improvement. 

A fully non-perturbative determination of the renormalisation factors of four-quark 
operators in the framework of the Schrodinger functional appears entirely feasible at this 
point, provided that one can overcome the well-known problem of the Eichten-Hill action, 
namely the exponential growth of statistical fluctuations at large Euclidean times [68]. 
Here the hope is that the methods described in refs. [27, 62] turn out to be as useful as in 
the simpler case of heavy- light bilinears. 

We have verified explicitly the expected mixing pattern in an extensive perturbative 
calculation at one loop. Thereby we have also obtained the NLO anomalous dimensions, 
which will be an important ingredient in future non-perturbative determinations of the 
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renormalisation factors. Furthermore, our perturbative calculation can be used to optimise 
the choice of renormalisation prescription in the forthcoming numerical simulations. 

As we have mentioned above, at the level of B^-B^ amplitudes the matching be- 
tween HQET and QCD requires to compute matrix elements of the operators Oyv+AA ^"^^ 
^SS+PP' ■^hich are mapped via tmQCD onto the operators CyA+AV ^^^ ^sp+pS' There- 
fore, as far as renormalisation is concerned, one is then faced with the task of computing the 
step-scaling functions for the relevant pair of operators, which renormalise multiplicatively, 

i.e. OyA+AV ^^*^ ^VA+AV + ^^SP+PS- 

The static approximation considered in this work only represents the lowest order 
of HQET, and hence all results for phenomenologically relevant quantities are subject to 
corrections in powers of the inverse heavy quark mass. While there are strategies in place 
which are designed for determining the leading 1/M corrections non-perturbatively [14], 
it is also possible to interpolate lattice results between the static approximation and the 
regime of relativistic quarks with masses around that of the charm quark. Our findings 
may serve to obtain high-precision results for B^ — B^ mixing amplitudes in the static 
approximation, which in turn are required to perform reliable interpolations to the physical 
6-quark mass. 
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A. Constraints from heavy quark spin symmetry and i^(3) spatial rota- 
tions on the mixing pattern 

We now describe the procedure followed to impose the constraints from heavy quark spin 
symmetry and cubic rotations. It applies identically to both the parity-even and the 
parity-odd sectors and we choose to describe it for the latter. It turns out that, in this 
particular case, there is no need for considering a maximal set of independent symmetry 
transformations, because the final constraints are already obtained by considering a finite 
spin transformation of the heavy fields, e.g. 

$1 : ii)h^ i>hl2ld. , i'h -^ V';i7273 , (A.l) 

and two lattice spatial rotations, e.g. 

<I>2 : 7^( 1^2) rotates the 1 axis onto the 2 axis, 

$3 : 1Z{ 2 — > 3 ) rotates the 2 axis onto the 3 axis, (A. 2) 



alone. The subspace spanned by (|2.q ) is not invariant under the set of transformations 
( |A.1| ). We hence give up temporarily Lorentz invariance (which, as we will see, will be 
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recovered naturally) and consider an enlarged basis containing eight operators, 



(A.3) 



which can generate, when properly combined, the original parity-odd basis ( |2.6D . The 
analysis performed in section |2| by using chiral symmetry can be carried over to ( [A.3| ), 
which can be accordingly shown to renormalise as 



O^n 



.± 



{l + 5^)0 



± 



(A.4) 



Here z^ are block diagonal matrices containing two (4 x 4) scale-dependent blocks while 
6 are block off-diagonal matrices, containing two (4 x 4) scale-independent blocks. The 
advantage of using ( A.3D is that the new basis is closed under ( [A.lj ) and ( [A.2| ), and the ma- 
trices $fc that implement the symmetry rotations can be constructed explicitly. Moreover, 
it should be observed that in order to preserve the renormalisation structure determined 
by chiral symmetry, the two matrices z and 6 have to satisfy the symmetry constraints 
independently, namely 



^± = $fcz±$ri , 



6^ = <^kS^<^^' 



The explicit form of the matrices ^k is easily found out to be: 



$1 



/O 1 

10 

1 

10 







\0 



o\ 



-10 

1 
1 0/ 



$9 



/l o\ 
00100000 
10 
10 
10 
10 
00000100 

\0 1/ 



$n 



(A.5) 



/l o\ 
10 
00010000 
00100000 
10 
00000100 
1 

\o 1 oy 



Once all the constraints are imposed one gets 



7. 



.± 



/ Zi Z2 Z2 Z2 

Z2 Zi Z2 Z2 

Z2 Z2 Zi Z2 

Z2 Z2 Z2 Zi 
















\ -4 



o\ 





-z' 



z. 

Z'2 Z\ 

z\ z\ 

z'o z[ J 



2 ^2 

1 ^2 ^2 

2 ^1 ^2 
2 ^2 
























-5'o -6'o -6'o 



6', 6[ 



S'2 
V ^2 



S'2 
^'2 



S'2 
S'2 



S'2 
5'2 
5'l 



-5i 62 §2 S2 \ 

-62 Si 62 S2 
-52 62 Si 62 
-62 62 S2 Si 



0/ 



"Notice that, due the constraints (2.3), the operators Osp+ps and Osp.ps are contained in the basis (A.3) 
^To simphfy the notation we neglect the superscript on 21, Z2, zj, Z2, 5i, 82, S'l, 82 despite the fact 
that these matrix elements are in general different in the S — ±1 sectors. 
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Now we go to a basis such that z is diagonal. A convenient choice is: 

{Qf,Qf + ^Qt,Rt,Rt,Qi + '^Qt,Qt-'iQt,Rt,Rt}, (A.6) 

where R- are some Lorentz non-invariant operators, the precise expression of which is 
irrelevant for the rest of the argument. By transforming into this basis it turns out that 
also the nontrivial blocks of 6^ are diagonal^. We have therefore arrived to the final 



conclusion that the basis ( 2.12 ) is closed under renormalisation, with the mixing pattern 
presented in ( 2.13| ). 



B. Integrals in the infinite volume theory 

The lattice contributions to the matching coefficients 'Vi-gfjat' commonly expressed in terms 
of a basic set of Feynman integrals in the momentum-space representation, have been 
calculated and cross-checked by independent authors [6,63-65]. In all cases, their evaluation 
has been pursued through Monte Carlo simulations (VEGAS), resulting in an average 
numerical precision of three digits. On the other hand, the matching constants -^gp i^f 
which provide the connection between the SF and the lattice scheme, have been calculated 



with a better accuracy, as explained in section 4.2. As a consequence, the uncertainty 
on the NLO coefficients of the SF anomalous dimensions is dominated by the lack of 
precision in the infinite lattice integrals. A possible way out would be running the Monte 
Carlo algorithms on faster computers and wait long enough for a couple of digits more. 
A more attractive alternative is to use an analytical trick to improve the quality of the 
results, obtaining at the same time some insight into the peculiar nature of the static lattice 
integrals. As an example, we consider 

3 1 



di 



\ r d'k 


H'-'%'^+4a!a. 


> sm -^ 

^^ 2 


4 

As = ^ sin^ k^ + i 


M=l 


fM=l 



^ 16 As 



Ai = ^sin^^, A2 = _^sin^fe^ + 4Af . (B.l) 

11=1 fj,=i 

The first term in the integral, which comes from the static propagator, diverges logarith- 
mically at A; = 0. This contribution is compensated by an opposite divergence of the 
subsequent terms, which brings the final result to a finite value, namely di ~ 5.46. In 
principle, di could be regularised through the usual lattice discretisation of the integration 

variables, 

27r A^ iV , , 

kfi -^ j^n^ , -— <n^< — . (B.2) 

If it were not for the 9(1 — A;^)-function, the integral would be expected to behave like an 
ordinary lattice integral, i.e. its convergence to the continuum would be determined by the 
asymptotic formula 

MN) = d, + ^ + |log(iV) + ^ + ^log(Ar) +0 (^) , (B.3) 



*The triplet {Q^ + 4Q2 1^1 5-^2 } corresponds to an eigenvalue with a three-dimensional associated 
subspace in the space of operators, and so does {Q3 — 2Qf, R^, Rf} - both for z and 5 . 



28 








1 1 1 1 1 


1 


1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 7: Convergence pattern oidi{N). 



where di represents the lattice version of di. The 0-function, which is non-zero inside a 
spherical domain, produces an explicit breaking of the hyper-cubic H{4:) symmetry, thus 
perturbing the convergence pattern. This effect can be better understood by defining 



/ = 



d^kOil-k^) — 
k^ 



i{N) = \ 



E 



Af2 



n 



n'* 



(B.4) 



and observing that the number of the lattice points that lie inside a 3-sphere S with radius 
i?s = N/2Tr increases irregularly when N ^ oo. The alternating excess or deficit of 
integration points gives rise to the oscillating behaviour characterising the dashed curve 
of Figure 4. In order to smooth di, we propose to regularise / in a way that formally 
restores the hyper-cubic symmetry. To this aim, we introduce an order parameter AVs, 
which provides a measure of the spherical symmetry breaking produced by the lattice 
discretisation, 

V^-V^iN) 



AFs(A^) 



Vs 



(B.5) 



Here Vs = A^'^/(327r^) is the volume of the above-mentioned 3-sphere, and Vs{N) represents 
the corresponding lattice volume, obtained by just counting the number of the lattice 
points that belong to the inner of S at fixed N. The recovery of spherical symmetry in 
the continuum limit implies that AVg vanishes when A^ — s- cxd. In addition, AVg describes 
a surface effect at large values of A^, and its rate of vanishing is therefore proportional to 
1/A^, up to fluctuations. We now define 

i'{N) = {l + aAVj:{N) + /3AVj:{Nf} i{N) , (B.6) 

where a and /3 are two real parameters to be suitably chosen. Although both / and I' 
diverge in the continuum limit, their difference vanishes. According to our considerations. 
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the fluctuations of AVs are expected to mimic the ones of /, and an appropriate and 
unique choice of a and /3 will provide a partial cancellation of the irregularities observed in 
the continuum approach of /', and consequently d\. The search of optimal values can be 
simply performed by hand, as far as only two parameters have to be tuned (in addition, /3 
multiplies a sub dominant contribution). We find, in particular, (a,/3) = (0.0838, —0.45). A 
plot of di, regularised according to this choice, is represented by the solid curve in Figure 4. 
A fit of the smoothed data against ( |B.3| ) allows to extract the value of d\ with much higher 
precision than the previous determinations. The procedure can be extended to all the other 
lattice integrals which contribute to the matching between the continuum and the lattice 
heavy-light operators. Their explicit expressions, reported in [6,29], will not be reproduced 
here, but a list of more accurate values, obtained with the method explained above, is given 
in Table || 



di 


Ji 


/ 


V 


c 


w 


5.4636(6) 


-4.8540(6) 


13.3503(6) 


-6.9230(24) 


4.5259(27) 


-1.20538(1) 



Table 2: Some basic infinite lattice integrals, which are needed to compute the matching coefficients 
"Vcf laf Their values have been determined according to the regularisation method explained in 
this Appendix. 



The finite lattice constants that contribute to the matching coefficients X 



(1) 



ref,lat 



are 



expressed in terms of those lattice integrals by a set of algebraic relations that have been 
first published in [6,29]. Their values are reported in Table ^. 



Bll 


-41.248(8) 


Dll 


2.5923(12) 


Dh 


-30.879(8) 


Drl 


2.4489(11) 


Dlr 


-37.843(6) 


Drr 


0.40179(3) 


dIr 


-37.843(6) 


D%, 


9.796(4) 


Drr 


-1.60717(1) 







Table 3: Values of the combinations of lattice integrals entering matching coefficients. 

The improvement integrals / , v and w , introduced in eq. ( [4.26| ), do not involve a 
factor of 6{1 — k"^). Hence, their numerical value can be computed with good precision 
through the blocking procedure described in [52]. We obtain: 

f = -3.6461(2) , 



w 



-6.7185(2) , 
0.82130(2) . 



(B.7) 



C. Tables. 

In this appendix we list our results for the coefficients r^g, as well as for the NLO anomalous 
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dimension, obtained for different discretisations and renormalisation schemes. We also take 
the opportunity to describe the blocking procedure [52] applied to determine the coefficients 
r^g (see also ref . [66] for another practical application in the context of the SF) . 

Here we apply this method at the level of the first two blocking steps, in order to 
eliminate the 0(a) cutoff effects in the data. Going beyond this level yields no benefit 
with double precision arithmetics, because of the numerical rounding that arises when 
subsequent cancellations of the signal are performed. Starting from the filtered data, we 
first check that the logarithmically divergent term in the one-loop renormalisation constant 
has the correct coefficient, which we indeed obtain with an average precision 

4,oh'k'^'^ = 1-000(1) , (c.i) 

cf. eq. (|]2l|). 

After having checked the form of the divergence, we remove it from the filtered data 
by subtracting explicitly a term 7^ '^ Mn(L/a). We then extract the finite part of the 
renormalisation constant r^Q according to the following procedure. The data are fitted 
with two different ansdtze, viz. 

n+l 



u=2 



Z(")(L/a)=A + ^(-j {5, + ain(L/a)} , n = l,2. (C.2) 



Rounding errors are modelled as suggested in [67]. We then fit in several intervals in L/a, 
always starting at the largest available value. Next we study the x^ P^r degree of freedom 
for each fit ansatz as a function of the number of values of L/a included, and find the 
minimum within the stability interval of the fit, thus obtaining two estimates [j^^glsp ^^^ 
['"fc'olsp ^°^ '^fc'o' ^^ ^^^^ t'^^olsp ^^ °^^ ^^^^ estimate. Then we consider a number of 
estimates for the uncertainty on r^Q, namely: the fit errors on [r-^glsp ^^'^ ['^fc^olsp' ^^^ 
difference [[j^^qIsp ~ I'^fc^olspl' ^^^ ^^^ fluctuation of A within the stability interval of each 
fit. The final uncertainty is taken to be the largest of all of them. 

Our final results for the finite coefficients of the renormalisation constants r^g are 



reported in Tables Q-^ below. In the determination of the NLO anomalous dimensions, we 
have derived their uncertainties by combining in quadrature the errors of r^Q and of the 
matching coefficients in eq. (4.25), cf. Table H. 
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action 



a 



1,0 



'2,0 



3,0 



'4,0 



-EH 


0.0 


1 


0.0 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


2 


0.0 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


3 


0.0 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


4 


0.0 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


5 


0.0 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


1 


0.5 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


2 


0.5 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


3 


0.5 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


4 


0.5 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 


-EH 


0.0 


5 


0.5 


-0.09022(12) 


-0.09029(4) 


-0.13884(12) 


-0.12270(4) 



Table 4: Finite parts of the renormalisation constants in the primed operator basis. Light quarks 
are regularised with the unimproved Wilson action (W); heavy quarks with the Eichten-Hill one 
(EH). Here Q — 0.0. Boundary conditions s are enumerated according to (3.4)-(3.8). 



action 



a 



1,0 



'2,0 



'3,0 



4,0 



SW-EH 


0.0 


1 


0.0 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


2 


0.0 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


3 


0.0 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


4 


0.0 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


5 


0.0 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


1 


0.5 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


2 


0.5 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


3 


0.5 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


4 


0.5 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 


SW-EH 


0.0 


5 


0.5 


-0.05219(12) 


-0.05005(13) 


-0.1046(13) 


-0.08498(15) 



Table 5: Finite parts of the renormalisation constants in the primed operator basis. Light quarks 
are regularised with the improved Wilson action (SW); heavy quarks with the Eichten-Hill one 
(EH). Here B — 0.0. Boundary conditions s are enumerated according to (3.4)-(|3~ 
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action 



a 



1,0 



2,0 



'3,0 



4,0 



-EH 


0.5 


1 


0.0 


-0.22302(8] 


-0.10054(7) 


-0.15083(5) 


-0.14615(10) 


-EH 


0.5 


2 


0.0 


-0.21290(8] 


-0.09671(7) 


-0.14767(5) 


-0.14187(10) 


-EH 


0.5 


3 


0.0 


-0.21290(8] 


-0.09566(7) 


-0.14767(5) 


-0.13734(10) 


-EH 


0.5 


4 


0.0 


-0.22302(8] 


-0.10054(7) 


-0.15083(5) 


-0.14615(10) 


-EH 


0.5 


5 


0.0 


-0.21290(8] 


-0.09356(7) 


-0.14767(5) 


-0.12828(10) 


-EH 


0.5 


1 


0.5 


-0.22723(8] 


-0.10475(5) 


-0.15504(10) 


-0.15036(6) 


-EH 


0.5 


2 


0.5 


-0.21711(8] 


-0.10092(5) 


-0.15188(10) 


-0.14608(6) 


-EH 


0.5 


3 


0.5 


-0.21711(8] 


-0.09987(5) 


-0.15188(10) 


-0.14155(6) 


-EH 


0.5 


4 


0.5 


-0.22723(8] 


-0.10475(5) 


-0.15504(10) 


-0.15036(6) 


-EH 


0.5 


5 


0.5 


-0.21711(8] 


-0.09777(5) 


-0.15188(10) 


-0.13249(6) 



Table 6: Finite parts of the renormalisation constants in the primed operator basis. Light quarks 
are regularised with the unimproved Wilson action (W); heavy quarks with the Eichten-Hill one 
(EH). Here Q = 0.5. Boundary conditions s are enumerated according to (3.4)-(3.8). 



action 



a 



1,0 



'2,0 



'3,0 



"4,0 



SW-EH 


0.5 


1 


0.0 


-0.1850(4] 


-0.06030(5] 


-0.11658(15) 


-0.10843(11) 


SW-EH 


0.5 


2 


0.0 


-0.1749(4] 


-0.05648(5] 


-0.11342(15) 


-0.10416(11) 


SW-EH 


0.5 


3 


0.0 


-0.1749(4] 


-0.05543(5] 


-0.11342(15) 


-0.09963(11) 


SW-EH 


0.5 


4 


0.0 


-0.1850(4] 


-0.06031(5] 


-0.11658(15) 


-0.10843(11) 


SW-EH 


0.5 


5 


0.0 


-0.1749(4] 


-0.05333(5] 


-0.11342(15) 


-0.09057(11) 


SW-EH 


0.5 


1 


0.5 


-0.1892(4] 


-0.06451(7] 


-0.12078(15) 


-0.11264(13) 


SW-EH 


0.5 


2 


0.5 


-0.1791(4] 


-0.06069(7] 


-0.11763(15) 


-0.10836(13) 


SW-EH 


0.5 


3 


0.5 


-0.1791(4] 


-0.05964(7] 


-0.11763(15) 


-0.10383(13) 


SW-EH 


0.5 


4 


0.5 


-0.1892(4] 


-0.06451(7] 


-0.12078(15) 


-0.11264(13) 


SW-EH 


0.5 


5 


0.5 


-0.1791(4] 


-0.05754(7] 


-0.11763(15) 


-0.09478(12) 



Table 7: Finite parts of the renormalisation constants in the primed operator basis. Light quarks 
are regularised with the improved Wilson action (SW); heavy quarks with the Eichten-Hill one 
(EH). Here B — 0.5. Boundary conditions s are enumerated according to (3.4)- (|3.8| ). 
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action 



a 



1,0 



2,0 



'3,0 



'4,0 



W-EH 


1.0 


1 


0.0 


-0.3596(10) 


-0.10601(14) 


-0.16560(14) 


-0.15650(23) 


W-EH 


1.0 


2 


0.0 


-0.3421(10) 


-0.10173(14) 


-0.16161(14) 


-0.15057(24) 


W-EH 


1.0 


3 


0.0 


-0.3421(10) 


-0.09975(14) 


-0.16161(14) 


-0.14540(24) 


W-EH 


1.0 


4 


0.0 


-0.3596(10) 


-0.10601(14) 


-0.16560(14) 


-0.15650(23) 


W-EH 


1.0 


5 


0.0 


-0.3421(10) 


-0.09579(14) 


-0.16161(14) 


-0.13509(24) 


W-EH 


1.0 


1 


0.5 


-0.3649(10) 


-0.11133(14) 


-0.17093(14) 


-0.16183(23) 


W-EH 


1.0 


2 


0.5 


-0.3474(10) 


-0.10706(14) 


-0.16693(14) 


-0.15589(23) 


W-EH 


1.0 


3 


0.5 


-0.3474(10) 


-0.10507(14) 


-0.16693(14) 


-0.15073(23) 


W-EH 


1.0 


4 


0.5 


-0.3649(10) 


-0.11133(14) 


-0.17092(14) 


-0.16183(23) 


W-EH 


1.0 


5 


0.5 


-0.3474(10) 


-0.10110(14) 


-0.16694(14) 


-0.14041(24) 



Table 8: Finite parts of the renormalisation constants in the primed operator basis. Light quarks 
are regularised with the unimproved Wilson action (W); heavy quarks with the Eichten-Hill one 
(EH). Here 9 ~ 1.0. Boundary conditions s are enumerated according to (3.4)-(3.8). 



action 



a 



1,0 



'2,0 



'3,0 



'4,0 



SW-EH 


1.0 


1 


0.0 


-0.3216(4) 


-0.06578(19) 


-0.13135(23) 


-0.1188(4) 


SW-EH 


1.0 


2 


0.0 


-0.3041(4) 


-0.06151(19) 


-0.12736(23) 


-0.1129(4) 


SW-EH 


1.0 


3 


0.0 


-0.3041(4) 


-0.05953(19) 


-0.12736(23) 


-0.1077(4) 


SW-EH 


1.0 


4 


0.0 


-0.3216(4) 


-0.06578(19) 


-0.13135(23) 


-0.1188(4) 


SW-EH 


1.0 


5 


0.0 


-0.3041(4) 


-0.05557(19) 


-0.12736(23) 


-0.0974(4) 


SW-EH 


1.0 


1 


0.5 


-0.3269(4) 


-0.07109(19) 


-0.13666(23) 


-0.1241(4) 


SW-EH 


1.0 


2 


0.5 


-0.3094(4) 


-0.06683(19) 


-0.13268(23) 


-0.1182(4) 


SW-EH 


1.0 


3 


0.5 


-0.3094(4) 


-0.06485(19) 


-0.13268(23) 


-0.1130(4) 


SW-EH 


1.0 


4 


0.5 


-0.3269(4) 


-0.07109(19) 


-0.13666(23) 


-0.1241(4) 


SW-EH 


1.0 


5 


0.5 


-0.3094(4) 


-0.06089(19) 


-0.13268(23) 


-0.1027(4) 



Table 9: Finite parts of the renormalisation constants in the primed operator basis. Light quarks 
are regularised with the improved Wilson action (SW); heavy quarks with the Eichten-Hill one 
(EH). Here 9 — 1.0. Boundary conditions s are enumerated according to (3.4)-(|3~ 
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6 


s 


a 


'1;SF / '1 


7S^"/Tr^'°' 


0.0 


1 


0.0 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


2 


0.0 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


3 


0.0 


-0.2087(4) + 0.02294(2)7Vf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


4 


0.0 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


5 


0.0 


-0.2087(4) + 0.02294(2)7Vf 


-0.5535(5) + 0.04876(3)iVf 


e 


s 


a 


'+;(!)/ '+;(o) 

73;SF /Ts 


'+;{!)/ '+;(o) 

74;SF /74 


0.0 


1 


0.0 


-0.0590(3) + 0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


2 


0.0 


-0.0590(3) + 0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


3 


0.0 


-0.0590(3) + 0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


4 


0.0 


-0.0590(3) +0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


5 


0.0 


-0.0590(3) +0.01371(2)iVf 


-0.1654(3) + 0.02478(2)iVf 



Table 10: The two loop anomalous dimensions of the diagonal basis in units of the corresponding 
universal one-loop coefficients. Here 9 = 0.0 and a = 0.0. 



e 


s 


Q 


'+;(!)/ '+;{o) 

7l;SF / '1 


7g^"/7r^'°' 


0.0 


1 


0.5 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


2 


0.5 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


3 


0.5 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


4 


0.5 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


0.0 


5 


0.5 


-0.2087(4) + 0.02294(2)iVf 


-0.5535(5) + 0.04876(3)iVf 


6 


s 


Q 


'+;(!)/ '+;{o) 

73;SF /73 


'+;{!)/ '+;(o) 

74;SF /74 


0.0 


1 


0.5 


-0.0590(3) + 0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


2 


0.5 


-0.0590(3) + 0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


3 


0.5 


-0.0590(3) + 0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


4 


0.5 


-0.0590(3) +0.01371(2)iVf 


-0.1654(3) +0.02478(2)iVf 


0.0 


5 


0.5 


-0.0590(3) +0.01371(2)iVf 


-0.1654(3) + 0.02478(2)iVf 



Table 11: The two loop anomalous dimensions of the diagonal basis in units of the corresponding 
universal one-loop coefficients. Here 6 — 0.0 and a — 0.5. 
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6 


s 


a 


'1;SF / '1 


I'A" /■/,*■■'" 


0.5 


1 


0.0 


0.1566(4) +0.00080(2)iVf 


-0.4690(5) + 0.04364(3)iVf 


0.5 


2 


0.0 


0.1287(4) +0.00249(2)iVf 


-0.5006(5) + 0.04555(3)iVf 


0.5 


3 


0.0 


0.1287(4) +0.00249(2)iVf 


-0.5092(5) + 0.04608(3)iVf 


0.5 


4 


0.0 


0.1565(4) +0.00080(2)iVf 


-0.4690(5) + 0.04364(3)iVf 


0.5 


5 


0.0 


0.1287(4) + 0.00249 (2) iVf 


-0.5265(5) + 0.04713(3)iVf 


e 


s 


a 


'+;(!)/ '+;(o) 

73;SF /Ts 


'+;{!)/ '+;(o) 

74;SF /74 


0.5 


1 


0.0 


-0.0327(3) + 0.01211(2)iVf 


-0.0364(3) + 0.01696(2)iVf 


0.5 


2 


0.0 


-0.0396(3) + 0.01254(2)A^f 


-0.0600(3) + 0.01839(2)iVf 


0.5 


3 


0.0 


-0.0396(3) + 0.01254(2)iVf 


-0.0849(3) + 0.01990(2)iVf 


0.5 


4 


0.0 


-0.0327(3) +0.01211(2)iVf 


-0.0364(3) + 0.01696(2)iVf 


0.5 


5 


0.0 


-0.0396(3) + 0.01254(2)7Vf 


-0.1347(3) + 0.02292(2)iVf 



Table 12: The two loop anomalous dimensions of the diagonal basis in units of the corresponding 
universal one-loop coefficients. Here 9 = 0.5 and a = 0.0. 



e 


s 


Q 


'+;(!)/ '+;{o) 

7l;SF / '1 


7s?'/7r^'°' 


0.5 


1 


0.5 


0.1681(3) +0.00010(2)iVf 


-0.4342(6) + 0.04153(3)iVf 


0.5 


2 


0.5 


0.1403(3) +0.00179(2)iVf 


-0.4658(6) + 0.04344(3)iVf 


0.5 


3 


0.5 


0.1403(3) +0.00179(2)iVf 


-0.4745(6) + 0.04397(3)iVf 


0.5 


4 


0.5 


0.1681(3) +0.00010(2)iVf 


-0.4342(6) + 0.04153(3)iVf 


0.5 


5 


0.5 


0.1403(3) + 0.00179(2)iVf 


-0.4918(6) + 0.04502(3)iVf 


e 


s 


Q 


'+;(!)/ '+;{o) 

73;SF /73 


'+;{!)/ '+;(o) 

74;SF /74 


0.5 


1 


0.5 


-0.0234(2) + 0.01155(2)iVf 


-0.0133(3) +0.01556(2)iVf 


0.5 


2 


0.5 


-0.0304(2) + 0.01197(2)iVf 


-0.0368(3) + 0.01699(2)iVf 


0.5 


3 


0.5 


-0.0304(2) + 0.01197(2)iVf 


-0.0617(3) +0.01850(2)iVf 


0.5 


4 


0.5 


-0.0234(2) +0.01155(2)iVf 


-0.0133(3) +0.01556(2)iVf 


0.5 


5 


0.5 


-0.0304(2) +0.01197(2)iVf 


-0.1116(3) + 0.02152(2)iVf 



Table 13: The two loop anomalous dimensions of the diagonal basis in units of the corresponding 
universal one-loop coefficients. Here 9 — 0.5 and a — 0.5. 
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e 


s 


a 


'+;{!)/ '+;(o) 

'1;SF / /I 


■^/i^''" 


1.0 


1 


0.0 


0.5321(3) - 0.02196(2)iVf 


-0.424(1) + 0.04090(8)7Vf 


1.0 


2 


0.0 


0.4840(3) - 0.01904(2)iVf 


-0.459(1) + 0.04304(8)iVf 


1.0 


3 


0.0 


0.4840(3) - 0.01904(2)iVf 


-0.475(1) + 0.04403(8)iVf 


1.0 


4 


0.0 


0.5321(3) - 0.02196(2)iVf 


-0.424(1) + 0.04090(8)iVf 


1.0 


5 


0.0 


0.4840(3) - 0.01904(2)iVf 


-0.508(1) + 0.04601(8)Aff 


e 


s 


a 


'+;{!)/ '+;(o) 

73;SF / '3 


'+;{!)/ '+;(o) 

74;SF /74 


1.0 


1 


0.0 


-0.0002(3) + 0.01014(2)iVf 


0.020(1) + 0.01351(8)iVf 


1.0 


2 


0.0 


-0.0090(3) + 0.01068(2)iVf 


-0.012(1) + 0.01549(8)iVf 


1.0 


3 


0.0 


-0.0090(3) + 0.01068(2)Aff 


-0.041(1) + 0.01721(8)Aff 


1.0 


4 


0.0 


-0.0002(3) +0.01014(2)Aff 


0.020(1) + 0.01351(8)Aff 


1.0 


5 


0.0 


-0.0090(3) + 0.01068(2)A^f 


-0.097(1) + 0.02065(8)Aff 



Table 14: The two loop anomalous dimensions of the diagonal basis in units of the corresponding 
universal one-loop coefficients. Here 9 = 1.0 and a = 0.0. 



e 


s 


Q 


'+;{!)/ '+;(o) 

7l;SF / /I 


7S?V7r^'°' 


1.0 


1 


0.5 


0.5467(3) - 0.02284(2)iVf 


-0.380(1) + 0.03824(8)iVf 


1.0 


2 


0.5 


0.4986(3) - 0.01993(2)iVf 


-0.415(1) + 0.04038(8)Aff 


1.0 


3 


0.5 


0.4986(3) - 0.01993(2)iVf 


-0.432(1) + 0.04137(8)Aff 


1.0 


4 


0.5 


0.5467(3) - 0.02284(2)iVf 


-0.380(1) + 0.03824(8)Aff 


1.0 


5 


0.5 


0.4986(3) - 0.01993(2)iVf 


-0.464(1) + 0.04336(8)Arf 


6 


s 


Q 


'+;{!)/ '+;(o) 

73;SF /73 


'+;{!)/ '+;(o) 

74;SF /74 


1.0 


1 


0.5 


0.0116(3) +0.00943(2)iVf 


0.050(1) +0.01174(8)iVf 


1.0 


2 


0.5 


0.0028(3) + 0.00997(2)iVf 


0.017(1) +0.01372(8)Arf 


1.0 


3 


0.5 


0.0028(3) + 0.00997(2)iVf 


-0.011(1) +0.01544(8)Arf 


1.0 


4 


0.5 


0.0115(3) + 0.00943(2)iVf 


0.050(1) +0.01174(8)A^f 


1.0 


5 


0.5 


0.0028(3) + 0.00997(2)iVf 


-0.068(1) + 0.01888(8)Aff 



Table 15: The two loop anomalous dimensions of the diagonal basis in units of the corresponding 
universal one-loop coefficients. Here 6 — 1.0 and a ~ 0.5. 
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